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Abstract. Measurements of elementary pp collisions are essential for under- 
standing heavy ion collisions. Results for pp collisions at 200 and 62.4 GeV are 
presented. At both energies NLO pQCD describes pion production well. The 
measured pion transverse single spin asymmetries are very large at 62.4 GeV 
O ' and are reasonably well described by models relying on pQCD at transverse 

momenta larger than 1 GeV/c. 
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One of the goals of the relativistic heavy ion program is to study the properties of 
matter at high temperature and high density. The explorations at the Relativistic 
Heavy Ion Collider (RHIC) indicate that that at cm. energies of 200 GeV per 
nucleon pair indeed such system is formed with novel properties as characterized 
through the particle production [ HI [21 [31 H] - Part of these conclusions rely on 
comparison to elementary pp collisions, where such nuclear effects should not be 
present. The RHIC experiments have also explored Au+Au and Cu+Cu collisions 
at 62.4 GeV, an intermediate energy between that at the SPS (18 GeV) and that 
where the bulk of the RHIC heavy ion program is performed at 200 GeV. The 62.4 
GeV was chosen to match the highest energy where data were obtained in the ISR 
experiments. Thus, to minimize the systematic uncertainty for the comparisons to 
the heavy ion data a brief run with (polarized) pp at 62.4 GeV was carried out 
to collect pp reference data by the BRAHMS, PHENIX and STAR experiments. 
Thus, at this point in time an extensive set of pp reference data have been collected 
at these two energies. 

In addition to the pp comparison data there is an extensive spin program at 
RHIC. It has long been envisioned that BRAHMS could contribute to these studies 
using the transversely polarized proton beams by measuring pions at large momenta 
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at forward angles. Part of run-5 and run-6 was dedicated to such measurements. 
There is current interest in understanding the large transverse single spin asymme- 
tries (SSA) observed at moderate to large values of xp in the energy range from 20 
to 200 GeV in the framework of different approaches to pQCD. The main theoret- 
ical efforts to account for the observed SSA have focused on the role of transverse 
momentum dependent (TMD) partonic effects in the structure of the initial trans- 
versely polarized nucleon[[5| and the fragmentation process of a polarized quark into 
hadrons[ [5]. Higher twist effects ("twist-3") arising from quark-gluon correlation 
effects have also been considered as a possible origin of SSA[[?1 |S]. To validate 
such an approach that is based on pQCD it is key to see how well pQCD describes 
inclusive data at high rapidities. 

This contribution presents results of inclusive pp reactions at 200 and 62.4 GeV, 
and comparisons to pQCD at both energies. Transverse spin results from 62.4 and 
200 GeV are also described and compared to model calculations. 

2. Inclusive pp measurements 

The data used for this analysis were collected with the BRAHMS detector sys- 
tem. The BRAHMS detector consists of two movable magnetic spectrometers, the 
Forward Spectrometer (FS) that can be rotated from 2.3° to 15° , and the Mid- 
Rapidity Spectrometer (MRS) that can be rotated from 34° to 90° degrees relative 
to the beam line, and several global detectors for measuring multiplicities and lu- 
minosity, and determining the interaction vertex, and providing a start time (TO) 
for time-of- flight measurement. 

The MRS is a single-dipole-magnet spectrometer with a solid angle of ~ 5msr 
and a magnetic bending power up to 1.2 Tm. Most of the data presented here 
were recorded at magnetic field settings of 0.4 and 0.6 Tm. The MRS contains two 
time projection chambers, TPM1 and TPM2 sitting in field free regions in front of 
and behind the dipole (D5). This assembly is followed by two highly segmented 
scintillator time-of-flight walls, one (TOFW) at 4.51 m and a second (TFW2) at 
either 5.58 m (90°setting) or 6.13 m (other angle settings). 

The FS consists of 4 dipole magnets Dl, D2, D3 and D4 with a bending power 
of up to 9.2 Tm. The spectrometer has 5 tracking stations Tl through T5, and par- 
ticle identifiing detectors: H2, a segmented time-of-flight wall, and a Ring Imaging 
Cherenkov Detector (RICH) [ [9] . Additional details on the BRAHMS experimental 
setup can be found in ref. [fTO], 

The minimum bias trigger used to normalize these measurements is defined with 
a set of Cherenkov radiators (CC) placed symmetrically with respect to the nominal 
interaction point and covering pseudo-rapidities that range in absolute value from 
3.26 to 5.25. This trigger required that at least one hit is detected in both sides of 
the array. 

The present analysis was done with charged particles that originated from col- 
lisions of polarized protons with interaction vertices in the range of ±40 cm. For 
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Fig. 1. Invariant cross section distributions for pion, kaons protons and anti-protons 
produced in p+p collisions at ^fs — 200 GeV at rapidities y = 2.95 (left panels) 
and y = 3.3 (right panels). In all panels, positive charged particles are shown with 
filled triangles and negative ones with open circles. The errors displayed in these 
plots are statistical. 



the 200 GeV data invariant cross sections were extracted in narrow (Ay = 0.1) 
rapidity bins centered at y = 2.95 and y = 3.3, respectively. Narrow rapidity bins 
are required to reduce the effects of rapidly changing cross sections in particular at 
higher pt- Each distribution is obtained from the merging of up to five magnetic 
field settings. The data are corrected for the spectrometer geometrical acceptance, 
multiple scattering, weak decays and absorption in the material along the path of 
the detected particles. Tracking and matching efficiencies for each of the 5 tracking 
stations in the spectrometer were calculated by constructing full tracks with only 
4 stations and evaluating the efficiency in the 5 th station. The overall efficiency is 
about 80-90% and is included in the extraction of the cross sections. Particles are 
identified by the RICH. The low momentum part of the proton spectra is measured 
using the RICH in veto mode. 

Data for y — 2.95 and y = 3.3 are presented for pions, kaons and protons in 
Figffl The pions exhibit a power law behavior, and at high pt the ratio ir~ /ir + is 
less than one indicating the increasing importance of the valence quarks. The p/p 
ratio is much smaller than 1 px indicating that in the case of fragmentation the 
gluons cannot dominate the particle production. 

The measured differential cross-sections are compared with NLO pQCD calcu- 
lations [[11] evaluated at equal factorization and renormalization scales, p = fip ~ 
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Fig. 2. Comparison of ir and K + invariant spectra at rapidity 2.95 to NLO pQCD 
calculations at 200 GeV. 



Hr = pr- These calculations use the CTEQ6 parton distribution functions [ [12] 
and a modified version of the "Knichl-Kramer-Potter" (KKP) set of fragmentation 
functions (FFs) [ |T5] referred to here as mKKP, as well as the "Kretzer" (K) set 
[ 1 14] . The KKP set includes functions that fragment into the sums ir + + n~ and 
K + + K~. Modifications were necessary to obtain functions producing the separate 
charges for both 7r and K. These modifications involve the following operations: to 
obtain the FFs producing 7r + , the functions that describe the fragmention of favored 
light quarks u, d into w° were multiplied by (1 + z) ( e.g. Dl = (1 + z)D n u with 
LP = ^D* +7r ) where z is the fraction of the parton momentum carried by the 

hadron, and D™ + the function fragmenting u quarks into positive pions. The func- 
tions fragmenting unfavored quarks u, d into ir° are multiplied by 1 — z. The same 
operation is done for 7r~, but this time the favored quarks are u and d. The FFs of 
strange quarks and gluons are left unmodified. Similar modifications were applied 
to obtain FFs into K + and K~ , but this time, the starting functions were the ones 
fragmenting u,U,s,s into the sum K + + K~ . Figure [5] shows that the agreement 
between the NLO calculations that include the mKKP FFs and the measured pion 
cross section is remarkable (within 20% above 1.5 GeV/c). Similar good agreement 
was obtained for neutral pions at y = [Q~5] and at y = 3.8 at RHIC. The 
agreement between the calculated and the measured kaon cross-sections is equally 
good. The difference between the mKKP and Kretzer parametrizations is driven by 
higher contributions from gluons fragmenting into pions. This difference has been 
identified as an indication that the gg and gq processes dominate the interactions at 
mid-rapidity [[15]. The present results indicate that such continues to be the case 
at high rapidity. The calculation that uses the Kretzer set underestimates the pion 
yields by a factor of ~ 2 at all values of pr while for positive kaons the agreement 
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Fig. 3. Invariant cross section for ir~ at rapidity 2.7 and 3.3 at 62.4 GeV. The 
curves are NLO pQCD calculations as described in the text. 



is good at low momentum but deteriorates at higher momenta. Additional details 
as well as comparisons to p and p spectra can be found in Ref. [ I17j . 

At 62.4 GeV , where the beam rapidity is 4.2, the spectrometer at forward 
angle samples produced particles that carries a significant fraction of the available 
momenta (31.2 GeV/c). Thus particle production is influenced by the kinematic 
limit. Data for identified charged hadrons were collected at 2. 3°, 3°, 4°, and 6°. 
Figure [3] shows differential cross sections for ir~ for rapidities 2.7 and 3.3. The cross 
sections changes rapidily with rapidity at high px where values up to 0.5 are 
probed. The data are compared to NLO pQCD in the same figure. The calculation 
are done in the same manner as for the 200 GeV using the KPP fragmentation 
function evaluated at /i = px scale. The calculation describes the overall magnitude 
and shape quite well, though at the highest rapidity there is a tendency for the 
calculation to fall below the data at the highest px- This may be in agreement 
with the analysis [ [18] of 53 GeV ir° data from the ISR at a fixed angle of 5° , 
comparable to the conditions for the present measurements (2.3°and 4°), but at 
larger xp where NLO pQCD results falls considerably below the data and with 
increasing discrepancy with xp. In contrast to the aforementioned paper we do, 
though, conclude that NLO pQCD gives a satisfactory description of the charged 
pion data at high rapidity. 



3. Transverse Single Spin Asymmetries 

The SSA is defined as a "left-right" asymmetry of produced particles from the 
hadronic scattering of transversely polarized protons off unpolarized protons. Ex- 
perimentally the asymmetry can be obtained by flipping the spins of polarized 
protons, and is customarily defined as analyzing power A^: 

1 (N+-CN-) 
N V{N+ + CN-y [> 
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where V is the polarization of the beam, C is the spin dependent relative luminosity 
(C = C+/C-) and is the number of detected particles with beam spin vector 

oriented up (down). The average polarization of the beam V as determined from 

BRAHMS Preliminary 




Fig. 4. An vs xjt for 7r + and tt~ at 62.4 GeV for positive and negative xp . Solid 
symbols are for ir + and open for ir~ . The curves are from the twist-3 calculations 
with (line) and without (broken) sea and anti-quark contribution. Predictions from 
the Sivers function calculation are shown as dotted lines. 

the on-line CNI measurements is about 60% for run 6 (62.4 GeV) and about 45% for 
run 5 (200 GeV). The systematic error on the An measurements is estimated to be 
20% including uncertainties from the beam polarization (~ 18%). The systematic 
error represents mainly scaling uncertainties on the values of An. 

The production rate of pions in given pp and xp bins are determined in a similar 
way as that for inclusive spectra, except that acceptance and effcicieny corrections 
are not performed since all of these will cancel out in the calculation of the quantities 
needed for determination of An- 

The analyzing power An for charged pions, An(tt + ) and An{tt~) at s/s = 62.4 
GeV as a function of xp are shown in Fig. Q] using data taken at 2.3° and 3° . 
The px range covered is from about 0.6 GeV/c at the lowest xp up to 1.5 GeV/c 
at the highest xp. The measured An values show strong dependence with xp 
reaching large asymmetries reaching up to ~ 40% at xp ~ 0.6. The An values are 
positive for tt + and negative for ir~ and with approximately the same magnitude. 
The asymmetries and their ^-dependence are qualitatively in agreement with the 
measurements from E704/FNAL. 

The analyzing power An for charged pions, An(tt + ) and An{k~) at y/s = 200 
GeV as a function of xp are shown in Fig. fusing data taken at 2.3° and 4° . The 
kinematic range is quite different from the 62.4 GeV data. The xj?-range covered 
is from 0.15 to 0.35, while the pt values covered by the measurements are much 
higher ranging from ~ 1 to 4 GeV/c. The left panel shows data from the 2.3° 
setting, while the right panel is for 4°. The latter for a given xp represents a higher 
Pt value than for the 2.3° data. Overall the asymmetries are smaller for the higher 
pr's. 
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Fig. 5. An vs xp for pions at 200 GcV. Solid symbols are for ir + and open for 
7r~measured at 2.3° (left panel) and 4° (right panel). The curves are from the twist- 
3 calculations with (line) and without (broken) sea and anti-quark contributions. 
Predictions from the Sivers function calculation are shown as dotted lines. 



Figure 0] compares An(tt) with a pQCD calculation in the range of pr > 
1 GeV/c using "extended" twist-3 parton distributions [ 7 J including the "non- 
derivative" contributions [ [TH1 [20] • In this framework, two calculations from the 
model are compared with the data: two valence densities (u v ,d v ) in the ansatz with 
and without sea- and anti-quark contribution in the model fit. The calculations 
describe the data within the uncertainties. 

For the calculations shown in the figure, the dominant contributions to SSAs 
are from valence quarks, while sea- and anti-quark contributions on SSAs are suffi- 
ciently small that the current measurements are not able to quantitatively constrain 
these contributions. The data are also compared with calculations using the Sivers 
mechanism which successfully describe FNAL/E704 An data. The calculations 
compared with the data use valence-like Sivers functions [ [21] [22] for u and d 
quarks with opposite sign. The fragmentation functions used are from the KKP 
parameterization [ [T3] , but the Kretzer fragmentation function [ Q3] gives similar 
results. The calculations shown with dotted lines in the figure underestimate An 
for both pt ranges, which indicates that TMD parton distributions are not sufficient 
to describe the SSA data at 62.4 GeV. The same kind of calculations are compared 
to the 200 GeV data in Fig O yielding similar conclusions as for 62.4 GeV 

4. Summary 

In summary, unbiased invariant cross sections of identified charged particles as func- 
tion of pt were measured at high rapidity in p+p collisions at y/s — 200 GeV and 
y^i = 62.4 GeV. NLO pQCD calculations reproduce reasonably well the produced 
particle (pions and kaons) distributions. The SSA for inclusive pions were mea- 
sured at forward rapidities at the same two energies. A twist-3 pQCD model of 
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An describes the xp and energy dependence reasonable well. There are theoretical 
challenges in describing the data at both lower energy and lower px- 
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